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Abstract 



We study the single-sink buy-at-bulk problem with an unknown cost function. We wish to 
' route flow from a set of demand nodes to a root node, where the cost of routing x total flow 

along an edge is proportional to f{x) for some concave, non-decreasing function / satisfying 
/(O) — 0. We present a simple, fast, combinatorial algorithm that takes a set of demands and 
constructs a single tree T such that for all / the cost f{T) is a 47.45-approximation of the 
optimal cost for that /. This is within a factor of 2.33 of the best approximation ratio currently 



c/3 . achievable when the tree can be optimized for a specific function. Trees achieving simultaneous 



O(l)-approximations for all concave functions were previously not known to exist regardless of 
computation time. 
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_ 1 Introduction 

Os 

(N 

^SJ ' Many natural network design settings exhibit some form of economies of scale that reduce the costs 

. when many flows are aggregated together. We may benefit from cheaper bandwidth when laying 

O I high capacity network links |AZ98| , reduced infrastructure costs or bulk discounts for shipping large 

amounts of goods together [SCRS97] , or summarization and compression of correlated information 
flows [KEW02] . These scenarios are known in the literature as buy-at-bulk problems. In a general 
buy-at-bulk problem we are given a graph and a set of demands for flow between nodes. The cost 
^ ' per unit length for routing a total of x flow along an edge is f{x) for some function /. To model 

. the economies of scale, we assume / is concave and monotone non-decreasing. 

We will focus on single-sink (or single-source) case, where all demands must be routed to 
a given root. When / is known, the problem becomes the well-studied single-sink buy-at-bulk 
(SSBaB) problem. SSBaB is A^P-hard — it generalizes the Steiner tree problem — but constant- 
factor approximations are known for any given / (e.g. [GMMOl] IGRlOj ). The special case where 
/ has the form /(x) = min{x,M} for some M (edges can be "rented" for linear cost or "bought" 
for a fixed cost) is known as the single-sink rent-or-buy (SSRoB) problem and has also received 
significant attention (e.g. [KMOQl EGRSIO] ) . 
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Buy-at-bulk algorithms produce trees that are heavily tailored to the particular function at 
hand, but in some scenarios / may be unknown or known to change over time. One setting where 
this arises is in aggregation of data in sensor networks. The degree of redundancy among different 
sensor measurements may be unknown, or the same network may be used for aggregating different 
types of information with different amounts of redundancy. In other situations rapid technological 
advancement may cause bandwidth costs to change drastically over time. Further, in the interest 
of simplifying the design process and building robust networks, it may be useful to decouple the 
problem of designing the network topology from that of determining the exact characteristics of 
the information or goods flowing through that network. In these settings it is desirable to find a 
single tree that is simultaneously good for all cost-functions, and from a theoretical perspective, 
the existence of such trees would reveal surprising structure in the problem. 

There are two natural objective functions which capture the idea of simultaneous approximation 
for multiple cost functions. Let T be the set of all concave, monotone non-decreasing cost functions 
satisfying /(O) = 0, /(T) be the cost of a routing tree T under function /, and be the optimal 
routing graph for /. Note that due to the concavity of / we may assume that is a tree. Let TZ 
be a randomized algorithm that returns a feasible routing tree T . First, we could try to minimize 
the quantity 

E7^[/(T)] 

sup ./^.N (1) 

which we call the oblivious approximation ratio. If the oblivious ratio is small, then TZ returns a 
distribution that works well in expectation for any /. However, there may be no sample from this 
distribution that works for everything: for any tree T there may be functions for which f{T) is 
expensive. 

To circumvent this problem we can work with the much stronger simultaneous approximation 
ratio. For a deterministic or randomized algorithm A that returns a tree T4 the simultaneous ratio 
of A is defined as: 



f{TA) 
sup -77;^ 



(2) 



A bound on the simultaneous ratio subsumes one on the oblivious ratio and proves there exists a 
single tree that is simultaneously good for all /. 

We emphasize that the distinction between the simultaneous and oblivious objectives is not a 
technicality in the objective but rather a fundamental difference and that the gap between these 
ratios can be large. Consider the problem of embedding arbitrary metrics into tree metrics, another 
problem that requires bounding the cost under many different functions (i.e. distortion of each edge). 
It is well-known that distributions over trees can achieve O(logn) expected distortion for all edges 
|FRT03] but that even for simple graphs like the n-cycle no single tree can do better than i}{n) 
distortion |RR98j . Therefore, the ratio between maximum expected distortion and the expected 
maximum distortion is r2(n/ log n) in this case. 

Goel and Estrin |GE03| introduced the problem of simultaneous SSBaB and gave an algorithm 
with an O(log-D) bound on the simultaneous ratio ([2|), where D is the total amount of demand. 
Goel and Post [ GP09j recently improved the oblivious ratio ([T]) to 0(1) for a large constant. Trees 
for which the simultaneous ratio was 0(1) were not known to exist regardless of computation time. 
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In this paper we give the first constant guarantee on the simultaneous ratio, resolving the major 
open question of Goel and Estrin and Goel and Post |GE03l[GP09j . Several aspects of our algorithm 
and analysis bear mentioning: 

• We achieve a simultaneous approximation ratio of 47.45. This is within a factor of 2.33 of 
the current best approximation for normal SSBaB of 20.42 ^GRlOj and substantially smaller 
than the 0(1) oblivious bound [GP09j . which we estimate to be around 15 million. 

• The algorithm is entirely combinatorial, in contrast to the result of Goel and Post [GP09] . 
which uses separation oracles to prove an 0(l)-oblivious distribution exists but reveals little 
of its structure. 

• Our analysis is short and simple, no more complex than the analysis of a normal SSBaB 
algorithm. 

• The runtime is only 0((t(n, m) + m + nlogn) log D) for a graph with n nodes, m edges, and 
D demand, where i(n, m) is the runtime of an SSRoB approximation. 

The algorithm is quite simple. We first find approximate trees for a set of rent or buy basis func- 
tions, prune this set to obtain a subset L of trees whose total rent costs are increasing geometrically 
while total buy costs are dropping geometrically, and then prove it suffices to approximate every 
tree in L. The set of bought nodes for each tree in L defines a series of tree layers, which we stitch 
together using light approximate shortest-path trees (LASTs) |KRY95j to approximate both the 
minimum spanning tree (MST) and shortest-path tree. Finally, we consider any layer in the tree. 
Using the geometrically changing costs and the properties of the LAST construction, we conclude 
that everything within the layer is an approximate MST, and everything outside approximates the 
shortest-path tree cost. 

1.1 Related Work 

The SSBaB problem was first posed by Salman et al. |SCRS97] . and the first general approximation 
algorithm was given by Awerbuch and Azar |AA97] . who used metric tree embeddings |Bar96] to 
achieve an O(log^n) ratio, later improved to O(logn) using better embeddings by Bartal |Bar98j 
and Fakcharoenphol et al. |FRT03] . Guha et al. |GMMOl| gave the first constant approximation, 
and follow-up work by Talwar [Tal02j . Gupta et al. |GKPR07] . Jothi and Raghavachari |JR04] . 
Grandoni and Italiano [Gl06], and Grandoni and RothvoB [ GRIP] has since reduced the constant 
to 20.42. Most recent algorithms for SSBaB (and several related problems) are based on the sample 
and augment framework of Gupta et al. [GKPROT] . Many algorithms using this framework have 
been derandomized by van Zuylen [vZ09] . 

The special case of SSRoB has also been extensively studied, often as a special case of the 
connected facility location problem. The first constant factor approximation was given by Ravi 
and Salman |RS99| as a special case of the traveling purchaser problem. Karger and Minkoff 
[KMOOj gave an alternate algorithm and introduced connected facility location. Gupta et al. im- 
proved the approximation to 9.01 [GKK+Ol] . Swamy and Kumar to 4.55 |SK04] . and Gupta et 
al. to 3.55 [GKPROT] . Gupta et al. [GST08] derandomized the 3.55-approximation to achieve a 
4.2-approximation. Eisenbrand et al. [EG RSlOj developed a randomized 2. 92- approximation, which 
recently improved to 2.8 using the 1.39-approximation for Steiner tree of Byrka et al. [BGRSIO] . 
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Since we employ the 2.8-approximation, and the claimed ratio does not currently appear else- 
where in the literature, we present the brief calculation deriving this value in the appendix. Both 
Williamson and van Zuylen [WvZ07| and Eisenbrand et al. [EGRSIO] independently derandomized 
the 2.92-approximation to achieve a deterministic 3.28-approximation. 

The problem of simultaneous approximation for multiple cost functions has been studied using 
both the oblivious and simultaneous objectives. Goel and Estrin |GE03] were the first to explic- 
itly pose the question of simultaneous approximations and gave an algorithm with an O(log-D) 
simultaneous guarantee. Prior to that Khuller et al. |KRY95j gave an algorithm to simultaneously 
approximate the two extreme cost functions f{x) = 1 and f{x) = x — a result which plays an im- 
portant role in this paper — and metric tree embeddings had been applied to SSBaB |AA97[ [FRT03] 
to achieve an O(logn) bound for the oblivious objective. Enachescu et al. |EGGM05] gave an 0(1) 
simultaneous guarantee for the special case of grid graphs with some spatial correlation. Goel and 
Post |GP09j proved that an oblivious guarantee of 0(1) is achievable for all graphs. Gupta et al. 
|GHR06j and Englert and Racke |ER09j have studied several generalizations of the problem where 
both the demands and function are unknown, and multiple sinks are allowed. In these settings the 
guarantee is generally O(logn) or 0(polylog n). 

2 Notation and Preliminaries 

Formally, we are given a graph G = (V, E) with edge lengths le for e G E', a root node r, and a set 
of demand nodes V <^V with integer demands d^. The total demand \s D = dy. We want to 
route dy flow from each w to r as cheaply as possible, where the cost of routing Xg flow along edge 
e is lef{xe) for some unknown, concave, monotone increasing function / satisfying /(O) = 0. Not 
knowing /, our objective is to find a feasible tree T minimizing supj f(T)/ f(TJ), where is the 
optimal graph for /. 

We first show that we can restrict our analysis to a smaller class of basis functions. Let 
e > be a small constant which will trade off the runtime and the approximation ratio, and 
K = [logi+,L>]. Fov < i < K, define Mi = (1 + ef, Ai{x) = min{a;,Mi}, and T* as the 
optimal tree for Ai. By the monotonicity and concavity of /, whenever Mi < x < Mj+i we have 
/(Mi) < f(x) < /(Mj+i) < (1 + e)f{Mi), so with a loss of only a factor of 1 -|- e we can interpolate 
between f{Mi) and /(Mj+i) and assume / is piecewise linear with breakpoints only at powers of 
1 -f e. A nondecreasing concave function that is linear between powers of 1 -|- e can be written as a 
nonnegative linear combination of {^i}o<i<_ft' by setting coefficients equal to the changes in slope: 
if the slope drops from 5i to at (1 -|- e)* it induces the term {6i — 5i+i)Ai{x). Now for a linear 
combination aiAi{x) and a tree T 

so it suffices to upper bound maxj ^j(T)/^j(T*). 

We now define some notation and subroutines that will be important for our algorithm. The 
problem of finding a good aggregation tree for the function Ai{x) = min{x,Mj} is an instance of 
the SSRoB problem, and we can find a A-approximate tree T,, where A is the best approximation 
ratio known, currently equal to 2.8 using the algorithm of Eisenbrand et al. [EGRSlO] and Byrka 
et al. [BGRSIO] . We will assume the algorithm is deterministic. If not (as in the case of the 2.8- 
approximation) we repeat it a polynomial number of times and pick the best tree, so we are close 
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to a A-approximation with very high probabihty. In this case, our simultaneous approximation 
algorithm will have some tiny probability of failure. 

The cost Ai(Ti) can be broken into two pieces, the rent cost and the buy cost, based on whether 
Ai is maxed out at Mj: 

Definition 2.1. For an aggregation tree Ti for cost function Ai with flow on edge e, the rent 
cost Ri and normalized buy cost Bi are defined as 

Ri = ^ ^ leAi{Xe) 

Note that edge costs composing Ri are weighted by the amount of flow they carry, but edges in 
Bi are not; they use unweighted edge costs. The total cost of Tj is given by ^i(Tj) = i?i + MiBi. 
The rent and buy costs also partition the nodes of Tj into two sets: 

Definition 2.2. The core Ci of tree Ti consists of r and all nodes spanned by bought edges and the 
periphery contains all vertices outside Ci. 

If we condition on the nodes in Ci then the rent-or-buy problem becomes easy: demands outside 
the core pay linear cost until they reach Q, so they should take the shortest path, whereas within 
Ci we pay a fixed cost per edge length, so the best strategy is to follow the min spanning tree. The 
cost Ri is therefore at least the sum of shortest path distances to Cj, while Bi is at least the weight 
of the MST of Q. 

In addition to the SSRoB approximation, we will also employ the light, approximate shortest- 
path tree algorithm of Khuller et al. [KRY95] : 

Definition 2.3 ( |KRY95] ) . For a > 1 and (3 > 1, an (a, /S)-light, approximate shortest-path or 
(a, /3)-LAST is a spanning tree T of C with root r such that 

• For each vertex v, the distance from v to r inT is at most a times the shortest path distance 
from V to r in G. 

• The edge weight of T is at most (3 times the weight of an MST of G. 

Khuller et al. show how to construct an (a, /3)-LAST for any a > 1 and /3 > Roughly, the 
algorithm performs a depth-first traversal of the MST of G starting from r, checking the stretch of 
the shortest path to each node. If the path to some v has blown up by at least an a factor, then 
it updates the tree to take the shortest path from v to r, adjusting other tree edges and distances 
accordingly. See the paper [ KRY95| for a full description and analysis. 

Finally, we define four parameters a, /3, 7, and 5 used by our algorithm whose values we will 
optimize at the end. 

• a > 1 is the approximation ratio for shortest paths used in our LAST. 

• /3 > is the corresponding approximation to the MST in the LAST. 

• 7 > 1 is the factor by which normalized buy costs Bi increase from layer to layer in our tree. 

• 5 > 1 is the factor by which rent costs Ri drop from layer to layer. 
We now turn to a more thorough explanation of tree layers. 
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3 Tree Layers 



In the normal SSBaB problem, the cost function is defined as f{x) = minj{(Tj + 5jx}, the cheapest 
of a collection of different "pipes" or "cables" given to the algorithm, each with an affine cost 
function aj + 6jX. It is common (e.g. [GMMOl] IGKPR07] ) to first prune these pipes to a smaller 
set with geometrically decreasing 5j's and geometrically increasing aj's and then build a solution 
in layers where each layer routes with a different pipe. 

We perform an analogous procedure. We would like to build our simultaneous tree T in a series 
of nested layers defined by the cores Cj of each tree Tj, so that the core of T under Ai is similar 
to Ci, but we have no guarantees on the relationships between different cores: Cj and Ck may be 
entirely disjoint except for r. However, we will show that as long as normalized buy costs Bi and 
are within a constant factor of each other, the same core can be used for both trees. Consequently, 
we are able to define nested layers by choosing one Ci for each order of magnitude of Bi. 

After finding A-approximate trees Tj for each Ai, we loop through the costs Bi and Ri, discarding 
i whenever Bi does not drop by 7 or Ri does not grow by 6. We are left with a subset L of the Tj 
where the Bi's are dropping by a factor of 7 and the iij's are growing by a factor of 6. The cores 
Cj for each i £ L will define the layers of our tree. Algorithm [1] describes the procedure in more 
detail. 

Algorithm 1: Finding tree layers 
Input: Graph G and demands T> 
Output: Set L and cores Cj for each i £ L 

1 for i ^ to do 

2 Tj A-approximate tree for Ai{x) 

3 for i ^ 1 to do 

4 if ^j(Tj_i) < ^j(Tj) then Ti ^ Tj_i 

5 for i ■(^ K — 1 down to do 

6 if ylj(Ti+i) < A,{Ti) then Tj ^ Tj+i 

7 for i ^ to do 

8 calculate Cj, Bi, Ri 

9 Lb ^ 

10 i? ^ 00 

11 for i ^ to do 

12 if Bi < ^B then 

13 Lb LbD {i} 

14 B ^ Bi 

15 T ^ 

16 i? ^ 00 

17 foreach i £ Lb in decreasing order do 

18 if Ri < ^R then 

19 L^LU{i} 

20 R -i^ Ri 

21 return L and Ci for each i £ L 



Intuitively, as it becomes more expensive to buy edges the optimum will buy fewer edges and 
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rent more. In the case of approximations, the progression becomes muddled because for some i 
the approximation guarantee may be tight while for i + 1 we may get lucky and find the optimum, 
resulting in both rent and normalized buy costs dropping. We first show that the monotonicity in 
buy and rent costs still holds as long as each Tj is better for Ai than both Tj_i and Tj+i. 

Lemma 3.1. After /me[3 of Algorithm{l\ for every i we have Bi > and Ri < iij+i. 

Proof. First we show that for each i, Ai{Ti) < min{^j(Tj_|_i), ^j(rj_i)}. After the loop on lines 
[SHU we have Ai{Ti) < Ai{Ti-i), and after the second loop on lines [SHSlwe have Ai{Ti) < Ai{Ti+i), 
so we only need to show that the second loop does not break the first condition. If the second loop 
updates Tj then Aj(Tj) will only shrink, and if it changes Tj-i it does this by setting Tj_i ^ Ti 
which preserves Ai{Ti) < Ai{Ti-i). 

Now consider AiiTf^) for any k. By definition Ai{x) < x and Ai{x) < Mj, so to upper bound 
Ai{Tk) we may assume edges within Ck pay Mj per unit length, which sums to MiBk, and edges 
outside Cfc pay linear cost, or i?^ total, implying Ai{Ti.) < Rj. + MiB).. Therefore 

R^ + MiB, = Ai{Ti) < Ai{T,+i) < Ri+i + MiBi+i 
=^ Mi{Bi - Bi+i) < Ri+i - Ri 

Similarly, 

Ri+i + Mi+iB,+i = A,+i{T,+i) < A+i{Ti) <Ri + M,+iBi 
=^ Ri+i - Ri < Mj+i(_Bj — -Bj+i) 

Combining the inequalities, 

MiiB, - B,+i) < R,+i -Ri< M,+iiBi - Bi+i) 

If Bi — -Bj+i < the inequality is false because Mj+i > Mi, so we conclude Bi > Bi^i. And using 
the first inequality, < Mi{Bi — -Bj+i) < Ri+i — Ri, so Ri < Ri+i. □ 

We need to show that we can restrict our attention to Tj for i £ L. Suppose i < k but Bi < ^B^. 
Using Lemma [3.11 observe that Ak(Ti) < Ri + M^Bi < Rk + SM^Bk < 5Ak{Tk). Note this is 
independent of the size of the intersection of the cores Ci and and any differences in routing. 
The following lemma generalizes this simple but key observation and proves that approximating 
each i G L is sufficient. 

Lemma 3.2. Suppose there exists a tree T and constants cb and cr such that for all i € L there 
exists a partition of the edges of T into two sets Tb^ and Tr. satisfying 

• ^0(7bJ < CBBi 

• AK{TR,)<CRRi 

then for allk £{0,..., K}, Ak{T) < max{cij7, Cij(5}AAfc(r^). 

Proof. Let k £ {0, . . . , K}. Let j = max{j G LbU < A;}. Either j = k or k was discarded due to j 
on lines II 1H14I because B^ > ^Bj, and either way Bj < jBj.. Now let i = min{i £ L\i > j}. Again 
i = j or j was pruned due to i on lines [T7VI201 and Ri < 6Rj . Applying Lemma 13.11 with i > j and 
j < k, we have Bi < Bj < ^Bj. and Ri < 6Rj < SR^. 
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This is sufficient to bound the cost of Ak{T): 

Ak{T) = Ak{TR^) + AkiTeJ < Ak{TrJ + MkAoiTsJ 

< CRRi + CBMkBi 

< CRdRk + CBlMkBk 

< max{cR6,CB^}Ak{Tk) < max{cR5,CBj}XAk{Tj^) 

The equahty follows because Tb, and Tr. partition the edges of T. The first inequality is because 
Ak{x) and MfcAo(x) both upper bound Ak{x), the second is by assumption, the third is from 
the derivation above, the fourth uses Ak{Tk) = i?^ + M^Bk, and the last is because T^, is a A- 
approximation. □ 

We will primarily assume that our SSRoB algorithm is a generic approximation, but Lemma 
13.21 can easily be improved to take advantage of an SSRoB algorithm with a stronger guarantee 
that separately bounds Ri and MiBi in terms of the optimal costs R* and MiB* . 

Corollary 3.3. Let T, cb, and cr he as in Lemma [Ql and suppose 

• Rr < f^RR* + f^BM.B* 

• MiBi < urR* + lyBM^B* 

then for all k, Ak{T) < max{ c^i^/Ur + crI'^R, crS^xb + CB'J'^BjAkiTj^). 
Proof. We change the inequalities in the proof above as follows: 

Ak{T) < CR6Rk + CB-fMkBk < CR6{fiRRl + fiR^hBl) + cbi{urRI + ubMuBD 

< max{ CR(5;UK + cr^vr, cr5^b + CB^i^B}Ak(T^) 

□ 

4 Constructing the Tree 

The construction of the tree itself is quite simple. We have a set of indices L and core sets Ci for 
i €z L. Starting with the largest i £ L, i.e. smallest Bi, and working downward, we connect each 
Ci to T, the tree so far, with a LAST. Algorithm [2] describes the procedure more formally. The 
notation G/T represents contracting T to a single node in G, and G[Ci] is the induced subgraph 
on Ci, so {G/T)[Ci] denotes first contracting T and then restricting to nodes in Cj. 

Lemma 4.1. The graph T constructed by Algorithmic is a tree and spans all demand nodes. 

Proof. Observe that G L, and Rq = 0, so Co covers all demands. Therefore after the last iteration 
T spans T). Each iteration only adds edges spanning new vertices, so no cycles are created. □ 

The tree may contain paths connecting Steiner nodes that carry no fiow. Such edges can be 
safely pruned or just ignored because they contribute nothing to the cost. 

All that remains is to define the partitions Tb^ and Tr. and prove the bounds needed in Lemma 
13.21 The set Tr. contains all edges present in T after connecting Cj, and Tr. contains the rest. Both 
cost bounds will follow easily from the geometrically changing costs: the cost Ao(TbJ is dominated 
by the cost of the Ci layer, an approximate MST, and Ak{Tr.) is dominated by the rent costs of 
the next layer, an approximate shortest-path tree. First, we bound the normalized buy cost of Tr.: 
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Algorithm 2: Constructing the tree 
Input: Graph G. Set L and accompanying Cj for each i G L 
Output: Aggregation tree T 

1 T ^ {r} 

2 foreach i €z L in decreasing order do 

3 T' ^ {a, /3)-LAST of {G/T)[Ci] with root T 

4 T ^TUT' 

5 return T 



Lemma 4.2. Lei i G L, and Tb^ be the tree T after the round when Ci is added. Then the edge 
cost j4o(TbJ is at most :^^Bi. 

Proof. The proof is by decreasing induction on i, i.e. in the order in which the layers are built. Let 
c be a constant to be chosen at the end. The base case is the largest i £ L, which is the smallest i 
such that Bi = 0. In this case, Ci = {r}, Tb^ = {r}, and the edge cost is 0. 

Now let i G L, = min{A; G L\k > i} be the previous (inner) layer, and suppose the edge cost of 
TB^. is at most cBk. By the construction of L, we know B^ < ^Bi, implying TB^. costs at most ^Bi. 
The cost of an MST of Ci in C is at most Bi, and Tg^ may already span part of Ci, so connecting 
the rest with an of {G /TB,.)[Ci\ costs at most Bi. Using an (a,/3)-LAST scales the cost by 

at most f3. 

The total cost of edges laid so far is at most ^Bi + /3Bi, so the proof is complete as long as 
cBi>^Bi + l3Bi. Setc=^: 

7 V 7/ 7-1 

□ 

Now we bound the rent costs of T/j. by a similar proof. 

Lemma 4.3. Let i £ L and Tr. = T/Tbh i.e. all edges outside ofTB^. Then the rent cost Ax(LijJ 
is at most -F-^^—rRi. 

Proof. We prove by increasing induction on i £ L (the reverse of Lemma 14. 2p that {Tr^ ) < cRi 
for some c to be determined. Since G L, and Tbq covers everything, the base case Trq costs too. 

For the inductive case, let i G L, A; = max{A: G L|/c < i} be the next (outer) layer, and Trj. 
have rent cost at most cRk. As before, note Rk < ^Ri, so Ak{Tbi,) < jRi- Tree Tg. spans Ci 
and possibly more, so if all demands outside Tg. took the shortest path from their sources to Tb^ 
the shortest-path cost would be at most Ri. However, the edges of Tr^, have moved some demands 
around, and by the time they reach the current layer they may be farther from Tg. than their 
original sources were. But by the triangle inequality the cost of sending all demands from their 
current locations to Tg^ via shortest paths is at most jRi + Ri, the cost of sending all flow in Tr^. 
back to its source and from there to Te. using shortest paths. The LAST algorithm guarantees 
a-approximate shortest paths, multiplying the cost by a. 

^We could allow Steiner nodes and use a Steiner tree approximation, but this would not improve the worst-case 
bound. 
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Consequently the total rent cost for Tr. is bounded by + a (f-Rj + Ri), which needs to be 



at most cRj. We can set c 



5-a-l ■ 



c ca f 1 a\ 5 — a — 1 a5 
c > - H — — + a <;=^ c 1 — - — r = c ; > a <;=^ c > 



6 6 \ 6 S J 5 - -5-a-l 

□ 

We note that Lemma [4.31 explains how we circumvent a major obstacle to an 0(l)-simultaneous 
approximation — the O(logn) distortion lower bound for embedding arbitrary metrics into tree 
metrics |Bar96| . If we needed to maintain distances between many pairs of nodes the task would 
be hopeless, but Lemma [4.31 shows that it suffices to preserve the distance of each node to the next 
layer, so the graph of distances to be maintained forms a tree. 

We can now complete the proof of our main theorem and choose the optimal parameters. 

Theorem 4.4. The tree T achieves a simultaneous approximation ratio o/ (1 + e)A(8 + 4\/5) using 
a X- approximation to SSRoB. In particular, 

• There is a randomized polynomial time algorithm that finds a 47-45 simultaneous approxima- 
tion with high probability. 

• There is a deterministic polynomial time algorithm that finds a 55. 58 simultaneous approxi- 
mation. 

• There exists a tree that is a 16.95 simultaneous approximation. 

Proof. Applying Lemma 13.21 with cb = (Lemma 14. 2p and cr = (Lemma 14. 3p . the final 

approximation ratio for an arbitrary cost function / is 

(1 + ejA max • 



7 — 1 6 — a — 1 

where the extra 1 + e comes from the approximation of / by a combination of ^j's. Now it is a 
simple matter of applying calculus to find the optimal values for a, 13, 7, and S. We set 

^^IWS /3 = 2 + V5 7 = 2 6 = 3 + V5 

for which = The derivation of these values is presented in the appendix. 

This gives us = 4(2 + i/S) so the simultaneous approximation ratio is 

(l + e)A(8 + 4^/5). 

Now, 

• Using the best randomized approximation A = 2.8, and the ratio is 47.45 with high probability. 

• Using the best deterministic approximation A = 3.28, and the ratio is 55.58. 

• If the algorithm is allowed to run in exponential time A = 1, and the ratio is 16.95. 
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□ 

The 2. 8- approximation of Eisenbrand et al. |EGRS10] actually provides a slightly stronger guar- 
antee on Ri and Bi, and we can use Corollary 13.31 to get a tiny improvement in the approximation 
ratio at the cost of a more complex derivation: 

Theorem 4.5. There is a randomized polynomial time algorithm that finds a 47.07 simultaneous 
approximation with high probability. 

Proof. Lemma 2 and Theorem 5 in [EGRSIO] prove that 

8n7 

E[R^] < 2R* + —^MiB* B[MiB,] < p{x + e)R* + pMiB* 

where p = 1.39 is the Steiner tree approximation ratio and x G (0,1] is a parameter. Applying 
Corollary 13.31 with 

.807 

PR = 2 pB = T^R = p[x + e) ^3= P 

the simultaneous ratio is bounded by 

(1 + e) max <^ + , — + 



Using 



5 — a — 1 7 — 1 ' x{5 — Q — 1) 7 — 1 



x = .5995 a = 1.5495 /3 = 7 = 2 5 = 2a + 2 

a — 1 



yields a simultaneous ratio of 47.07. The parameters are derived in the appendix. □ 

We leave as an open question the problem of exploiting Corollary 13.31 to substantially improve 
the ratio. 



4.1 Runtime 

Let m) be the running time of our SSRoB approximation on a graph with n vertices and m 
edges, which must be at least Q.{n) to write down the output. When e is constant, running the 
SSRoB approximation for each i takes 0{t{n,m) log D). Subsequent loops in Algorithm [1] take 
0(n log D) = 0{t{n, m) log D). 

For each of the O(log-D) iterations of Algorithm [2] we need to do a graph contraction and run 
the LAST algorithm, which requires computing the MST and shortest path trees. The computation 
of the shortest path tree takes 0{m + nlogn) and dominates the other steps. Combining the two 
algorithms, the total time is 0{{t{n,m) + m + nlogn) log D). 
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5 Open Problems 



We have answered the open questions posed by Goel and Estrin and Goel and Post |GE031 [GP09] . 
but there are several avenues for further work. Our simultaneous ratio of 47.45 already surpasses 
many algorithms for normal SSBaB and is only a factor of 2.33 away from the best. It would be nice 
to eliminate this gap or, alternately, prove that a gap exists between the approximation achievable 
for fixed / and the best simultaneous ratio. We know of no lower bounds on what simultaneous 
ratio may be possible, so any progress in this direction would also be interesting. Generalizing the 
settings in which 0(1) simultaneous ratios are possible would be interesting, but may be unlikely 
given that one must contend with lower bounds for metric tree embedding |Bar96 j and multi-sink 
buy-at-bulk [And04j. 
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A Derivation of optimal parameters 

Here we derive the optimal values for the parameters used in the proofs of Theorems 14.41 and 14.51 
Parameters for Theorem \4-4\ We need to minimize the expression 

r 1 

max < , > [6) 



7 — I'J — a — 1 
The easiest parameters to fix are 7 and 5. For 7: 



d_ 

d^ 



7-IJ (7-1) 
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For 6: 

a6^ 



d_ 

d5 



5-a-l 



25{5 - a - 1) - 
" {5-a-lY =° 



^ 52 - 2a(5 - 25 = 5(5 - 2a - 2) = ^ (5 = 2a + 2 
Plugging 7 = 2 and 5 = 2a + 2 into ([3]), = A/3, and 

a52 a(2a + 2)2 



5-a-l (2a + 2) - a - 1 



4a(a + 1) 



so ([3]) is now max{4/3, 4a(a + 1)}. The constraints on a and /3 require /3 > |KRY95j . so one 
term blows up if the other shrinks. To minimize the maximum set the two expressions to be equal: 

a + 1 2 ^ ^ 

p = = a(a + 1) =^ a(a — 1) = 1 =^ a — a — 1 = =^ a = — - — 

a — 1 2 

Using a = , we get 

^ a + 1 3 + V5 (3 + ^/5)(l + ^) 8 + 4^5 ^ , /- 
P = 7 = ^ = ^ = -. = 2 + V5 

and 5 = 2a + 2 = 3 + \/5. □ 
Parameters for Theorem \4.5\ We need to minimize 

r 2a52 />(x + e)/37' -807052 ^ 2 . 
[0 — a — 1 7 — 1 x(5 — a — 1) 7 — IJ 

For 7 and 5, both expressions inside the max function are minimized exactly as above in Theorem 
14.41 with 7 = 2 and 5 = 2a + 2. Also as in Theorem 14.41 (3 can be set to the minimum allowed 
by the constraints. Plugging in these values and simplifying, expression [J] becomes 

f p(x + e) .807a p 
4(a + 1) max <^ 2a + ^-^ + ^ 



a — 1 X a — 1 

Set the two expressions inside the maximization to be equal, and solve the resulting quadratic 
equation in x to get (for e = 0): 



_1 a(a-l) /a2(a-l)2 .193a(a - 1) 1 

^ ~ 2 p y p'^ P 4 



using that x > 0. 

The problem is now to minimize 



4(a + l)(2a + ^)=4a(a + l) + M^ + 4(a + l)Ja2-:l^+ 



a — 1 J a — 1 y a — l(a — 1)2 

This expression is unwieldy to optimize analytically but when p = 1.39 it achieves a minimum of 
about 47.068 for a ^ 1.5495. When a = 1.5495, x ^ .5995. □ 
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B Approximation ratio for SSRoB 



The current-best algorithm for SSRoB [EGRSIO] depends on the approximation ratio for Steiner 
tree, which has recently been reduced to 1.39 [BGRSIO] . The improved SSRoB ratio does not 
currently appear in the literature, so we include the calculation for completeness. See the original 
paper for details. 

Theorem B.l ( [EGRSlOi [BGRSIO^ . There is a 2. 8- approximation for SSRoB. 

Proof. By the proof of Theorem 6 in [EGRSlO] . there is an SSRoB algorithm with expected cost 

piMB* + {x + e)R*) + 2R* + 0.807 

X 

where p is the approximation ratio for Steiner tree and x is a parameter in (0, 1]. Set the coefficients 
of R* and MB* to be equal and solve the resulting quadratic equation in x. For p = 1.39, choosing 
X = .5735 gives an approximation ratio of 2.80. □ 
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